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The expansion and gravitational collapse of 
spherical balls of dust are investi- 
gated. Three cases are considered: positive, zero and 

negative curvature of the space within the cloud of 
dust. In each of these three cases the interior geometry 
is matched to an exterior Schwar zschild geometry. Also, 
in each case, the mass m seen by an observer at infinity 

is found to be positive definite and related to the den- 

3 

sity p ar.d radius r of the dust via m = (4iT/3)pr . The 
radius and density of the dust are functions of time. 

The method developed in this paper can be applied to 
thick shells of dust surrounding a collapsing or sta- 
tionary body, giving simple models of supernovae. The 
envelope expands leaving the central body as a remnant. 
These results are compared to those of Newtonian physics. 
Some effects of rotation are discussed. 
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I. 


INTRODUCTION 


In Newtonian mechanics one investigates the expansion or 
collapse of a homogenous spherical ball of dust by integrat- 
ing Newton's equations. The same problem in general rela- 
tivity is more complicated; it usually involves solving 
Einstein's field equations in two coordinate patches and 
then fitting the geometries together. The patching problem 
is often the most difficult part of the problem. However, 
in the case of expanding or collapsing dust, an interesting 

theorem due to Beckedorff and Misner^ and Lindquist and 
2 

Wh' ;ler reduces the patching problem almost to the sim- 

3 

piie.'.ty of Newtonian mechanics. This theorem states that: 
a necessary and sufficient condition for two 4-geometries 
to join smoothly (at the interfact) is that particles at 
the interface follow simultaneously the geodesic laws of 
motion for the two separate geometries. The method is 
applied in ref. 3 for the case of positive curvature of 
the space within the dust. Through the use of this theorem 
one can match the geometries without matching the metrics. 

In this paper we consider the cases of zero and 
negative curvature as well as (for completeness) positive 
curvatuie. One obje* tive is to write the solution in such 
a way that it can be used in a subsequent paper concerning 


2 


the gravitational collapse of rotating balls of dust. 

In section II, Newton's equations for the gravitational 
collapse of a ball of dust are written down and integrated. 

This is doen to facilitate the interpretation of the rela- 
tivistic resuits for the cases of zero and negative curva- 
ture. For the case of positive curvature, it is often more 

convenient to interpret the solution via the initial value 

4 

equations at the point of time symmetry. In sections III 
and IV geodesics are investigated. The geometries are 
matched and the results are discussed physically in sections V and 
VI. In section VII, the effect of rotation is considered. 

II. NEWTONIAN COLLAPSE 

To facilitate the interpretation of later results, a 
short Newtonian investigation of the expansion or collapse 
of a dust ball is given here. Newton's equaaions for the 
outermost dust particle of the dust ball are 

-2 

Mr’ = -Mm r . (1) 

Here M denotes the geometrized mass^ of the dust particle 
m, the mass of the entire dust ball r, the radius of the 
ball, and the dot denotes differentiation with respect to 
time. Integration of Eq. (1) yields 

f 2 = -K + 2m r ^ (2) 

n 


3 


whore -K is integration constant equal to twice the total 
n 

energy of the particle per unit mass. 

If the cases of positive, 2 ero and negative are 
investigated separtely, the solution takes the simple form 


. 2 

r " r o ain ( r /2) 

t = (h) (r^ /2m) ^ ( r| -sin r,) 


(3a) 


for 


0 < K = 2m/r (def. of r ) , 
no o 


r = 


Bt 


2/3 


B 


3 


S n/2 for K = 0, 
n 


(3b) 


and 


for 


r = r sinh (r/2) 
o 

t = ’h) (r /2m) (sinh r\ - r;) (3c) 

o 


0 > K = -2m/r . 
n o 


For K >0 the dust expands and then recontracts. In 
n 

the other two cases the dust only expands or contracts. In 
general relativity, the cases of positive, zero and 
negative correspond to positive, zero and negative curvature 


- 4 - 


of the space within the dust. Also the expressions (3) 
become identical with those of general relativity if the 
time is replaced by the proper time. 


III. INTERIOR SOLUTION 

The space interior to a homogeneous distribution of 
dust is often described by the Robnrtson-Walker^ Metric: 


2 2 2 r 21-2 r 2 22 2 2 

ds = -dT + a (T) 1 + (k/4) U I dU + U d8 + U sin Od 

l J L 

(4) 

Here K =1, 0, -1 corresponds to positive, zero negative 

curvature of the space within the dust. 

The quantity a, which is proportional to the luminosiL 

radius of the dust ball, must satisfy the field equation 

oo oo 

G = 8nT . 


In a synchronous frame 6 on a three dimensional spacelike 
surface, this equation takes the form 


(3) 


R + K 2 = 16 np 


where is the scalar curvature of the spacelike surface, 

K is the extrinsic curvature related to the second funda- 
2 


mental form K. . via 
13 


K. 


i n 2 i j 

- K J K. . , 
13 


■2 * ("\ 

and p is the energy density of the dust. 
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Direct calculation yields 


f ^ R = 3k a ^ and = 3 a 

Here a denotes differentiation with respect to T. 
Consequently the above field equation takes the form 

k a ^ = (16 tt p/3) - (a/a) ^ 

5 

The Robertson-Wa Iker Metric (4) is often used to 
describe a homogeneous and isotropic universe (a Friedman 
universe) . The average density p necessary to yield a 
given Hubble constant H = a/a depends on the curvature 
constant k. A larger density is necessary for positive 
curvature than for zero or negative curvature as can be 
seen from inspection of the above equation. 

When these three cases are treated separately, we 
obtain results very similar to those of Eq. (3) : 

2 

a = a sin ( r /2) 
o 

T = (a /2) ( r\ -sin r) 
o 

3 

a = (8tt/3) pa 
o 

for k = 1, 


a 

3 

a 


aT 


2/3 


6 npa 

for k = 0, 


(5a) 


(5b) 


6 


0 


= a^ sinh* - ( ^/2 ) 

T = (a /2) (sinh - - (5c) 

o 

o = (3r/3) t .ci 3 
o 

for k = -1, 

Here p is the time dependent density of the dust while 

a and a ore mteqration constants, 
o 

In the space described by this metric there are geodesics 
for which U = constant and T is the proper time along the 
geodesic. This is not surprising since the dust was assumed 
to be co-moving in the deviation of this metric. 


IV. EXTERIOR SOLUTION 

The space exterior to a sphere of dust can be described 
by the Schwa rzsc'ni Id metric: 


J 2 2 2 2 2 22 2.2..: 

ds = -A at + B dr + r dS + r sm * dc 


where 



-2 

B 


1-2 m 


-1 

r 


Here the constant m is the mass seen by an observer at in- 
finity. The connection between this macs and the density 
of the dust can be obtained by matching the interior and 
exterior solutions (patching) . 
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An especially simple way to obtain the geodesic equa- 
tions for radial motion of a dust particle is this geometry 
(6) is via the variotional principle 

o=6 ds (7) 

These geodesic equations arc 

- 2 g ' 

A t' 1 =0 (3) 

L -d 

and 

_ A~ 2 r ' ~ , + (A 2 ) (t') 2 "| + A -4 (r') 2 = 0 (9) 

L r i_ _i 

where the prime denotes differentiation with respect to 

the proper time 

2 2 2 - 

T ( dT*" = -ds ) and A =1 -2m r 1. Integration of Eq. (6) 
yields 

A 2 t' = K (10) 

where K is an integration constant. From Eq. (9) , t' can 
be eliminated via substitution of Eq. (10) yielding 

’a" 2 (r') 2V + K 2 A -4 (A 2 ) ' = 0 (11) 

I J 

after some manipulation. 

Integration of Eq. (11) yields 

(r') 2 = K 2 + K A 2 = K 2 + K - 2m K r’ 1 , (12) 


6 



where K is i. n in Leg rot ion cons ton t. For large r and small 
velocities, Eq. (2) should reduce to o Newtonian expression 
(2). The two expressions (12) and (2) agree in this limit 
if K - -1. Hence, Eq. (12) 1 i ke Eq . (2) contains only one 

adjustable constant K. The first integral (12) can also be 
obtained directly from the line element by dividing all 
terms by the proper time and eliminating t' via Eq. (10). 
This latter method gives K = -1 automatically. 

The solution of the general relativistic equations of 
motion takes a simple form when the three cases (l-K**) 
positive, zero and negative are treated seaprately. The 
solution is very similar to that for the Newtonian case: 


2 

r = r sin (-/2) 
o 

T = w) (r 3 /2m) ‘ ( rj -sm **) 

o 

2 -1 

for 0 < 1 - K =2mr 

o 


r = 


T 


2/3 


2 

for K 


= 1 , 



and 


r = r sinh^ i' r /2) 
o 

T = ( l i) (r^ /2m) 2 (sinh r, - r) 

2 -1 
for 0 > 1 - X = -2m r 


(13a) 


(13b) 


(13c) 
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The expressions (3) and (13) differ only in that the time 
in Eq. (3) is replaced by the proper time in Eq. (13). 

The solution (13) describes the motion of a test 
particle in gravitational field of a spherically symmetric 
body ni mass m. When tho maBcivo body ic rompoRPd of 
gravitationally collapsing dust, there is a connection 
between the mass seen by an observer at infinity and the 
density of the dust. This connection is obtained by 
matching the interior and exterior solutions. 

V. MATCHING OF GEOMETRIES AND DISCUSSION 

Application of the theorem discussed in section I, 
reduces the problem of matching four-geometries of that 
of matching eodesics. In section III and IV the geodesics 
can be matched by inspection . The motion of the same par- 
ticle is describable in two different coordinate patches. 

Thus, in each description, the proper time of a particle 
at the inteface between the two geometries and the proper 
circumference of the ball of dust should agree. 

In view of the -above discussion, comparision of Eqs. 

(5) and (13) yields 

Tl + (k/4) U 2 T 1 a U = r 
L o J o o o 

o is 

a = (r j /2m) 2 (14a) 

o o 

for k = 1 or -1 
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and 

au = ? 

o 

, (14b 

for 

k = 0. 

Thus, if the conditions (14) are satisfied, the four-geome- 
tries will join smoothly at the interface. These conditions 

Q 

(14) when combined with Eqs. (5) and (13) yield 

m = ( 4 tt/ 3 ) p r 3 (15) 

for each case. Hence, the general relativistic expression 
for the mass of a 'Homogeneous ball of dust is identical with 
the well known Newtonian expression in terms of the density p 
and radius r of the dust ball. Also, in each case the proper 
time required for collapse from a finite radius to the gravi- 
tational radius is finite. 

To an observer at infinity, the radius of the dust 
ball always appears to be greater than the Schwarzschild 
radius. Thus, an observer is unable to see the later stages 
of the collapse because of the time dilation. Since no 
energy escapes from the dust ball, the mass seen by this 
observer is constant. On the other hand, if the ball of 
dust (with r > 2m) is initially expanding it will recon- 
tract only for the case k = 1. 


11 


Using the method described above there con be constructed 
o universe which has sections with negative as well os posi- 
tive curvature and with all of these sections fitting together 
smoothly. A simple example can be constructed as follows: In 

the above calculations we have obtained the geometry exterior 
and interior to a collapsing dust ball with negative curva- 
ture within it. If a positive curvature Friedman universe is 
then patched to the large r portion of the Schwarzschild re- 
gion, we have a universe. The method used to patch these two 
regions together is identical with that used in this paper 
to patch the dust ball to the Schwarzschild region. Ic is 
not known whether the resulting universe is open or closed. 

As on aid to visualizing the situation, consider a 
closed universe at the instant of time symmetry containing 
a Schwarzschild region with a positive curvature Friedman 
region at each end. Such a universe is shown in Fig. 1. 

The universe with a negative curvature section is very 
similar, it has a Schwarzschild region with a positive curva- 
ture Friedman region on one end and a negative curavature 

Friedman region on the other. Such a universe can expand in 
one region and contract in another. 

Other results of physical interest can ha obtained via 
differentiation of Eq. (13) with respect to the proper time, 
or substitution of Eq. (13) into Eq. (12), yielding 



(16a) 


r ' = (2m/r) (1 - r r^ 1 ) 

for k = 1 

r = ( 2m/r) 

for k = 0 , 
and 

k -1 \ 

r ' = ( 2m/r) (1 + r r Q ) 

for k = -1. 


.2 4 2 

the line element (6) r = A (r') 


A 2 ♦ (r') 2 ^ 1 . 


This velocity r ' is related to the velocity r = dr/dt via 

Thus , 

the velocity seen by an observer at infinity (who remains 
a constant distance from the center of the dust ball) is 


r = 


2 r -1 

h i 2mr 


for k =1, 



r = A (2m 


r" 1 ) 


for k = 0 
and 


2 r —1 — 1 ~ 1^5 P — 1 “|k 

r = A ' 2m r + 2m r i 1 + 2m r ! 

L o j L o J 

for k = -1. 

In Newtonian mechanics, there are three cases: (1) 

dust expands and then recontracts, (2) the dust expands 


(16b) 


(16c) 


a t 


(17a) 


(17b) 


(17c) 


the 


13 


continuously with the expansion velocity r approaching zero 

as r approached infinity, (3) the dust expands continuously 

but the expansion velocity approaches a finite value as a 

r approaches infinity. Inspection of Eq. (17) shows that in 

. . 9 

general relativity these three cases correspond to positive , 
zero, and negative curvature respectively.. In the latter 

— 1 H —1 ^2 

case r approaches (2m r ) (1 + r Q 2m ) as r approaches 

infinity. 

When combined with Eq. (17) , Eq. (10) has a simple physi- 
cal interpretation. For k = 1 and r = r Q , Eq. (10) takes the 
form 

-1 h 

dT = (1 - 2m r ) 2 dt (18) 

o 

which is the well known relation between time intervals of 
the Schwarzschild solution. For k = 0 and k = -1, the dust 
expands to infinite radius. Consequently, the expression 
(10) can be compared with the results of special relativity 
since the Schwarzschild geometry is asymptomatically flat. 
According to special relativity, the time interval dT 
measured by an observer sitting on a moving particle is 
shorter than the time interval dT seem by a stationary 
observer i.e. 

dT = (1 - r 2 )^ dt. (19) 
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Substitution of the value of r for r approaching infinity 
yields 


t' = (dt/dt) = 1 

for k = 0 , 

and 

t = ( 1 + 2m r ■*" ) ^ 
o 

for k = -1 where t = dt/dt. 

These expressions are identical with Eq. (10) for 
r -*■ infinity and k = 0 and -1 respectively. Thus, in 
certain regions the results of Eq. (10) can be obtained 
by simple physical arguments. 

VI. DISCUSSION OF ASTROPHYSICAL APPLICATIONS 

The results of this paper have astrophysical ap - 
plications . These results allow a general relativistic 
description of astrophysical systems which expand and 
never reccntract as well as those which recontract e.g. 
exploding stars ( super novae ) , galaxies, or clusters of 
galaxies at the stage where the pressure and rotation 
are negligible. The negative curvature case corresponds 
to an unbound system in Newtonian mechanics, one which 
never recontracts if it is initially expanding. The 
zero cur mature case corresponds to the transition be- 
tween bound and unbound systems (the zero total energy 
case) ; it just barely expands to infinity using up all 


( 20 ) 


( 21 ) 


its kinetic energy in the process. Although it is possible, 
it is very unlikely in the k = 0, -1 case that the velocity 
vectors will be arranged such that the system will collapse. 
Consequently, in a physical situation, the k = 0, -1 cases 
are more likely to be expanding than contracting. 

Using the method discussed in this paper, one can 
construct general relativistic supernovae models. The 
expanding envelope of the supernovae corresponds to the 
expanding negative curvature solution given here. For 
example, to construct such a model of a supernovae, (col- 
lapsing or static) to an exterior Schwarzschild solution, 
then patch this to an expanding negative curvature solu- 
tion, and finally patch this to another exterior 
Schwarzschild solution. For such a model, the envelope 
blows off, leaving a remnant behind. The mass of this 
remnant m r is given by 

m r = (4u/3) pr 1 3 

where p is the density and r^ is the inner radius of the 
envelope. In this way, we can treat the expansion or 
collapse of thick shells of dust. 
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VII. EFFECT OF ROTATION 

An important question in astrophysics is: Does 

rotation stop collapse or does collapse crush rota- 
tion 11 ? In order to study this question, consider a 
test particle on the surface of the collapsing ball of 
dust. In Newtonian mechanics one learns that if the 
particle has sufficient tangential velocity, it will 
remain in orbit around the collapsing dust. If the 
particle remains in orbit, it has a point of closest 
approach to the center of the collapsing dust ball. 

At this point the outward radial acceleration is equal 
to or greater than zero. The former case corresponds 
to circular motion. 

The general relativistic equations of motion of 
the test particle can be obtained from the variational 
principle (7) yielding 


2 ' ' 

(r <p ) = 0 

2 ' ' 

(A^ t ) = 0 


2B 2 r" + (B 2 ) r (r') = 2r ((£*) 2 - (A 2 ) r (t') 2 . 


(22a) 

(22b) 

(22c) 
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Integration of Eqs. (22a) and (22b) 

r 2 cp ' = i (23a) 

A 2 t' = K (23b) 

and substitution into Eq. (22c) yields 

2 2 2 2-1 - ° 2 -4 

2B r" + (B ) (r') = 2^ r - 2mr K A . (23c) 

At the test particle's point of closest approach to 
the center of the dust ball, the radial velocity r' vanishes 
and the r.-.dial acceleration r'' is equal or greater than 
zero. This implies that 

2 2-4 

l * m r K A (24) 

On the other hand, the world line of the particle must re- 
main within the light cone. This is true only if 

2 2 2 -2 

K < K r A * (25) 

at the point of closest approach. These two conditions (24) 

and (25) are consistent only for r 5 3m. 

Thus, for r s 3m, rotation of the test particle about 

the center of the dust ball can stop the collapse of the 

1Z 

particle . The particles can remain in orbit abound the 
collapsing dust ball. However, for 3m > r > 2m the gravi- 
tational attraction overpowers the centrifugal force. 
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Thus, in this region the rotation cannot stop the collapse. 
This is because the expression for centrifugal force devi- 
ated from the Newtonian one since the particle world line 
must remain within the light cone; the local velocity of 
the particle cannot exceed that of light. 

A collapsing star can shed some of its angular momen- 
tum via the above mechanism and collapse toward its 
Schwarzschild radius leaving planets in orbit about it. 

But for 3m > r > 2m the planets cannot remain in orbit. 

The centrifugal force is overpowered by the gravitational 
attraction and everything collapses. 
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Figure Caption 


Figure 1: Closed inhomogeneous universe at the moment of 

of time symmetry. This universe is obtained 
by patching positive curvature Friedman regions 
to both enas of a Schwarzscmia region. 



